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' We compute the one-loop supersymmetric (SUSY) contributions to the weak 

CN ■ 

^ . charges of the electron (Qh/)) proton {Q^r), and cesium nucleus (Q^y) in the Minimal 

Supersymmetric Standard Model (MSSM). Such contributions can generate several 
percent corrections to the corresponding Standard Model values. The magnitudes 
. of the SUSY loop corrections to Q^^ and Qy^ are correlated over nearly all of the 

, MSSM parameter space and result in an increase in the magnitudes of these weak 

O ■ P 

' charges. In contrast, the effects on Q)^ are considerably smaller and are equally 

o ■ 

■ likely to increase or decrease its magnitude. Allowing for R-parity violation can lead 
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to opposite sign relative shifts in Qy^ and Q^^, normalized to the corresponding 
Standard Model values. A comparison of Qy^ and Qyy measurements could help 



^ ' distinguish between different SUSY scenarios. 
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I. INTRODUCTION 

The search for physics beyond the Standard Model (SM) of electroweak and strong 
interactions is a primary objective for particle and nuclear physics. Historically, parity- 
violating (PV) interactions have played an important role in elucidating the structure of 
the electroweak interaction. In the 1970's, PV deep inelastic scattering (DIS) measure- 
ments performed at the Stanford Linear Accelerator Center (SLAC) confirmed the SM 
prediction for the structure of weak neutral current interactions [1]. These results were 
consistent with a value for the weak mixing angle given by sin^ Ow 1/4, implying a tiny 
\^(electron) Xy4(quark) neutral current interaction. Subsequent PV measurements - per- 
formed at both very low scales using atoms as well as at the Z-pole in e^e~ annihilation - 
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have been remarkably consistent with the results of the SLAC DIS measurement [1]. 

More recently, the results of cesium atomic parity- violation (APV) [2] and deep inelastic 
V- {p-) nucleus scattering [3] have been interpreted as determinations of the scale-dependence 
of sin^^^vK- The SM predicts how this quantity should depend on the momentum transfer 
squared (g^) of a given process ^ . The cesium APV result appears to be consistent with the 
SM prediction for 0, whereas the neutrino DIS measurement implies a +3a deviation 
at Iq^l ~ 20 (GeV/c)^. If conventional hadron structure effects are ultimately unable to 
account for the NuTeV "anomaly" , the results of this precision measurement would point 
to new physics. 

In light of this situation, two new measurements involving polarized electron scattering 
have taken on added interest: PV Moller (ee) scattering at SLAC[4] and elastic, PV ep 
scattering at the Jefferson Lab (JLab)[5]. In the absence of new physics, both measurements 
could be used to determine sin^ 9^ at the same scale ^ 0.03 (GeV/c)^] - falling between 
the scales relevant to the APV and neutrino DIS measurements - with comparable precision 
in each case ^ (Asin^ ~ 7 x 10~^). Any significant deviation from the SM prediction for 
sin^ 9w at this scale would provide striking evidence for new physics, particularly if both 
measurements report a deviation. On the other hand, agreement would imply that the most 
likely explanation for the neutrino DIS result involves hadron structure effects within the 
SM. 

In this paper, we analyze the prospective implications of the parity-violating electron 
scattering (PVES) measurements for supersymmetry (SUSY). Although no supersymmetric 
particle has yet been discovered, there exists strong theoretical motivation for beheving that 
SUSY is a component of the "new" Standard Model. For example, the existence of low- 
energy SUSY is a prediction of many string theories; it offers a solution to the hierarchy 
problem; and it results in coupling unification close to the Planck scale. In addition, if R- 
parity is conserved (see below), SUSY provides an excellent candidate for cold dark matter, 
the hghtest neutrahno (see Ref. [7] for a review). The existence of such dark matter is 
required in most cosmological models. In light of such arguments, it is clearly of interest to 

^ The weak mixing angle and its evolution are renormalization scheme-dependent. Here, we use the MS 
scheme for the SM and the DR scheme for supersymmetric extensions of the SM. 

^ In practice, the PV ep experiment will actually provide a value for sin^ 0w{<l^ = 0), as discussed in Ref. 
[6]. 
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determine what insight about SUSY-if any-the new PVES measurements might provide. 

In the simplest version of SUSY - the Minimal Supersymmetric Standard Model (MSSM) 
with conserved R-parity [8] - low-energy precision observables experience SUSY only via tiny 
loop effects involving virtual supersymmetric particles. The requirement of baryon minus 
lepton number (B — L) conservation leads to conservation of the R-parity quantum number, 
Pr — (— 1)2'S'+3(B-L)^ where S denotes spin. Every SM particle has Pr — +1 while the 
corresponding superpartner, whose spin differs by 1/2 unit, has Pr — —1. Conservation 
of Pr implies that every vertex has an even number of superpartners. Consequently, for 
processes like ee — > ee and ep — > ep, all superpartners must live in loops, which generate 
corrections - relative to the SM amplitude - of order (a/Tr) (M/M)^ ~ 10~^ (where M 
denotes a SM particle mass and M is a superpartner mass). Generally speaking, then, low- 
energy experiments must probe an observable with a precision of few tenths of a percent 
or better in order to discern SUSY loop effects. Low-energy charged current experiments 
have already reached such levels of precision, and the corresponding implications of these 
experiments for the MSSM have been discussed elsewhere[9]. 

In the case of PV ee and elastic ep scattering, the precision needed to probe SUSY loop 
effects is roughly an order of magnitude less stringent, owing to a fortuitous suppression of 
the SM PV asymmetries, A^^j. At leading order in q^, the A{e) x V{f) contributions to 
Aj^ji are governed by Q{y, the "weak charge" of the target fermion, /. The weak charge of 
a particle / is defined as the strength of the effective A{e) x V{f) interaction: 

G 

^EFF = -^Q^eT^Tse/T/./ ■ (1) 

At tree-level in the SM the weak charges of both the electron and the proton are suppressed: 
Q^r — —QiY — 1~4 sin^ 6\y 0.1. One-loop SM electroweak radiative corrections further re- 
duce this tiny number, leading to the predictions Q^^r — —0.0449 [6, 10] and Q^r = 0.07 16 [6]. 
The factor of > 10 suppression of these couplings in the SM renders them more transparent 
to the possible effects of new physics. Consequently, experimental precision of order a few 
percent, rather than a few tenths of a percent, is needed to probe SUSY loop corrections. 
(Theoretical uncertainties associated with QCD corrections to are considerably smaller 
[6, 10]). 

In analyzing these SUSY loop contributions to Qy^, we carry out a model-independent 
treatment, avoiding the choice of a specific mechanism for SUSY-breaking mediation. While 



4 



most analyses of precision electroweak observables have been performed using one or more 
widely-used models for SUSY-breaking mediation, the generic features of the superpartner 
spectrum implied by such models may not be consistent with precision data [9]. Conse- 
quently, we wish to determine the possible impact of SUSY on the two PVES measurements 
for all phenomenologically acceptable choices of the MSSM parameters, even if such choices 
lie outside the purview of standard SUSY-breaking models. In doing so, we follow the spirit 
of Ref. [11], where a similar analysis of SUSY loop effects in v (p)-nucleus scattering was 
performed. 

In the case of PV electron scattering, we find that the magnitude of SUSY loop effects 
could be as large as the proposed experimental uncertainties for the Q\y and Q^, measure- 
ments (8% and 4%, respectively [4, 5]). Moreover, the relative sign of the effect (compared to 
the SM prediction) in both cases is correlated - and positive - over nearly all available SUSY 
parameter space. To our knowledge, such correlation is specific to the MSSM (with R-parity 
conservation), making it a potential low-energy signature of this new physics scenario. We 
also find that the SUSY loop effects on Q^y , the weak charge of the cesium atom measured 
in APV, is much less pronounced. Thus, the present agreement between the experimental 
value for and the SM prediction does not preclude the presence of relatively large effects 
in the PV electron scattering asymmetries. 

We also investigate a scenario where Pr is not conserved. We find that, in contrast to 
the Pij-conserving SUSY, the relative sign of the effect (compared to the SM prediction) is 
always negative for and can have either sign for Q^, with positive sign being somewhat 
more likely than the negative sign. The potential magnitude of the effects are considerably 
larger than those generated by SUSY loops. In principle, then, a comparison of Q^r and 
can potentially establish whether or not R-parity is violated within a SUSY extension of the 
SM. Having an answer to this question would have consequences reaching beyond the realm 
of accelerator physics. For instance, if Pr is violated in PVES, then lepton number is not 
conserved, thereby implying that neutrinos have Major ana masses and making neutrinoless 
double beta decay possible (see e.g. Ref. [12]). R-parity violation also renders the lightest 
super symmetric particle unstable, thus eliminating SUSY dark matter, which has significant 
implications for cosmology [7]. 

Our discussion of these points is organized as follows. After briefly reviewing the Minimal 
Supersymmetric Standard Model in Section II, we discuss the structure of the one-loop 
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radiative corrections to in Section III and the tree-level Pr- violating contributions in 
Section IV. The analysis of the prospective implications of the parity-violating electron 
scattering measurements for supersymmetry is presented in Section V. We conclude in 
Section VI. Appendix A lists the counterterms for the effective PVES Lagrangian of Eq. (1) 
necessary for renormalization of the one-loop radiative corrections to the weak charges. In 
Appendix B we explicitly prove that gluino loops do not contribute to Q^. In Appendix 
C we give complete expressions for all process-dependent one-loop SUSY corrections to the 
ep and ee scattering; expressions for process-independent contributions are given in the 
appendices of Ref. [11]. 

II. MSSM PARAMETERS 

The content of the MSSM has been described in detail elsewhere[8], so we review only a 
few features here. The particle spectrum consists of the SM particles and the corresponding 
superpartners: spin-0 sfermions (/, which include sneutrinos z/, charged sleptons Z, and up- 
and down-type squarks u and d), spin-1/2 gluinos (^), spin-1/2 mixtures of neutral Higgsinos 
(^°_2), the bino (5), and the neutral wino (IV^), collectively called neutralinos (Xi-4); 
and spin-1/2 mixtures of charged Higgsinos [H^) and charged winos (VT^), collectively 
called charginos (xi^2)- addition, the Higgs sector of the MSSM contains two doublets 
(up- and down- types, which give mass to the up- and down-type fermions, respectively), 
whose vacuum expectations and are parameterized in terms oi v — \jv\ + v\ and 
tan/3 = Vu/vd- Together with the SU(2)i and U(l)y couphngs g and g' respectively, v 
is determined from a, M^, and G^, the Fermi constant extracted from the muon lifetime, 
while tan (3 remains a free parameter. The MSSM also introduces a coupling between the two 
Higgs doublets characterized by the dimensionful parameter ji. The complete set Feynman 
rules for the MSSM, which take into account SUSY breaking and particle mixing, are given 
in Ref. [13]. 

Degeneracy between SM particles and their superpartners is lifted by the SUSY-breaking 
Lagrangian, which depends in general on 105 additional parameters. These include the 

SUSY-breaking Higgs mass parameters; the electroweak gaugino masses Mi^2\ the gluino 
mass Mg\ the left- (right-) handed sfermion mass parameters M| (M| ); and left- right mix- 
ing terms M? which mix Jl and Jr into mass eigenstates /i,2. In our analysis, we take 
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the sfermion mass matrices to be diagonal in flavor space to avoid large flavor-changing 
neutral currents. We also set all CP-violating phases to zero. One expects the magnitude 
of the SUSY-breaking parameters to he somewhere between the weak scale and ~ 1 TeV. 
Significantly larger values can reintroduce the hierarchy problem. 

Theoretical models for SUSY-breaking mediation provide relations among this large set 
of soft SUSY-breaking parameters, generally resulting in only a few independent parameters 
at the SUSY-breaking or GUT scale[14]. Evolution of the soft parameters down to the 
weak scale introduces flavor- and species-dependence into the superpartner spectrum due to 
the presence of Yukawa and gauge couplings in the renormalization group (RG) equations. 
According to the model-independent analysis of Ref. [9], however, generic features of this 
spectrum implied by typical SUSY-breaking models and RG evolution may conflict with the 
combined constraints of low-energy charged current data, M^, and the muon anomalous 
magnetic moment unless one allows for non-conservation of Pr. In light of this situation, we 
adopt here a similar model-independent approach and do not impose any speciflc relations 
among SUSY-breaking parameters. To our knowledge, no other model-independent analysis 
of MSSM corrections to PV observables has appeared in the literature, nor have the complete 
set of corrections to low-energy PV observables been computed previously (see, e.g. Ref. [15] 
for a study within minimal supergravity and gauge mediated models of SUSY breaking). 



With higher-order corrections included, the weak charge of a fermion / can be written as 



where and Qf are, respectively, the weak isospin and the electric charge of the fermion 
/. The quantities ppv and Kpv are universal in that they do not depend on the fermion 
/ under consideration. The correction Aj, on the other hand, does depend on the fermion 
species. At tree-level, one has ppv — 1 — kpv and A/ = 0, while at one-loop order, these 
parameters are 



III. RADIATIVE CORRECTIONS TO Q 



Qw = Ppv 2/3 - AQfKpv sin^ 9w + A/ , 



(2) 



Kpv 



Ppv 



1 + <5^SM ^ ^^SUSY ^ 
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FIG. 1: Types of radiative corrections to parity- violating electron scattering: (a) Z boson self- 
energy, (b) Z — 7 mixing, (c) electron anapole moment contributions, (d) vertex corrections, and 
(e) box graphs. External leg corrections are not explicitly shown. 

where the SUSY contributions to ppy, Kpy, and A/ are denoted in the above equation by 
the corresponding superscript. In general, the corrections (5p, 5k, etc. depend on g^, and 
in particular, the g^-dependence of Kpy defines the scale-dependence of the weak mixing 
angle: sin^ = Kpy(q'^) sin^6'vy, with sin^ being evaluated at some reference scale 

(usually = M|). 

The precise definitions of sin^ /«py(g^), etc. depend on one's choice of renormalization 
scheme. We evaluate the SUSY contributions using the modified dimensional reduction 
renormalization scheme (DR) [16] and denote all quantities evaluated in this scheme by a 
hat. In DR, all momenta are extended to d = 4 — 2e dimensions, while the Dirac algebra 
remains four- dimensional as required by SUSY invariance. The relevant classes of Feynman 
diagrams are shown in Fig. 1. Note that all gauge boson self-energies contribute only to ppy 
and Kpy while all non-universal box diagrams as well as vertex and external leg corrections 
are combined in A/. The counterterms for the effective PVES Lagrangian in Eq. (1) in the 
DR scheme are given in Appendix A. 

The Z boson self-energy contribution (Fig. la) simply rescales the leading order am- 
plitude. Its effect is naturally combined with the counterterm 5(5^ from Eq. (A2) into 
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Ppv- 

5G tizm 1 nw(o) ^ 

where the Z boson self-energy is evaluated at = 0. This is an appropriate approximation 
in this case because the momentum transfer in ee and ep scattering will be much smaller 
than the masses of the particles that appear in the loop graphs. The error is of the order 
|g^|/M| ~ 10~^, which is negligible. The quantity 5y^ denotes the sum of electroweak 

vertex, external leg, and box graph corrections to the muon decay amplitude, which must 
be subtracted when the neutral current (NC) amplitudes are normalized to G^. 
The graphs Fig. lb,c contribute to upv in Eq. (2). The expression is: 

«py = 1 + - 2 + 4c'Fl + — , 5 
s 

where F^[q'^) is the parity-violating electron-photon form factor, which - at = - is 
known as the anapole moment of the electron (see Eq. (12) below) ^. It should be noted 
that in the MSSM one has n|^^'^(g^) ~ g^, so there is no singularity at g^ = in the above 
equation [17] . The SM contribution contains a singularity that is canceled by a corresponding 
singularity in the anapole moment contribution (Fig. Ic). Since in the following we consider 
only the new physics contributions, this issue is irrelevant (a complete treatment of the SM 
contributions is given in Ref. [6, 10, 18]). The shift Ss"^ in = 1 — = sin^^iy(M|) arises 
from its definition in terms of a, G^, and Mz'- 



-2-2 

s c — 



V2G^MUl - Af) ' 

Af - n' I o^^n^7(o) t^zzjMl) n^^(o) ^ 

where ni^^(g^) = n^^/g^. Writing Af = Af^^^ + Ar^^^^ one has 

"■^SUSY ^ C A -SUSY (j\ 

In computing the SUSY corrections to the weak charges one must decide which value for 
to use. Since ^Sgusy has already been absorbed into Kpv one must determine from 



^ Note that our definition of Kpv, Eq. (5), includes the anapole form-factor of the electron F|(g^), which 
may be absent in other definitions appearing in literature (see e.g. Ref. [6]). 
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Eq. (6) using the SM radiative corrections only. The corresponding value extracted using 
only a, G^, and Mz is [12]: 

= 0.23120 ±0.00018 . (8) 

In order to incorporate constraints from existing precision data (see Section V), it is 
useful to introduce the oblique parameters S, T, and U [19]: 

\ New 



2a2 



-^Re|n..(o)-n..(Ml) + £!^ 



T 



1 



2s^ 



n^^(M|)-n^^(o)J+n^^(M|) 

\ New 



nzz(o) + — nz^(o) - Tiww{^) 



jj _ 4£^/ nw(o) - nw(M^) ., nzz(M|)-nzz(o) 

- ^ Ml 



a 



+ 2cs 



nz^(M|) - nz^(o) , ,2^,,{Mf 



New 



(9) 



where the superscript "New" indicates that only the new physics contributions to the self- 
energies are included. Contributions to gauge-boson self energies can be expressed entirely 
in terms of the oblique parameters 5", T, and U in the hmit that Mnew ^ Mz- However, 
since present collider limits allow for fairly light superpartners, we do not work in this limit. 

Consequently, the corrections arising from the photon self-energy (H-y^) and 7-Z mixing 
tensor (H^^) contain a residual g^-dependence not embodied by the oblique parameters. 
Expressing ppy and Kpv in terms of S',T, and U we obtain: 



a 



Ml 



SUSY 



+ 



+ 



a 



SUSY 



;2 77e/„2\SUSY 



(10) 



where Ad is the SUSY contribution to the difference between the fine structure constant 



SUSY 



As 



and the electromagnetic coupling renormalized at = Mz'- Act = [a{Mz) — ol 
noted above, we take — > in our analysis. 

The non-universal contribution to the weak charge is determined by the sum of the 

renormalized vertex corrections Vy^ in Fig. Id [see Eq. (A6)] and the box graphs 5^^^ in 
Fig. le [see Eq. (C15)] : 

>^f-9iVX + g\v4 + 5t^, (11) 
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where 



Qyj^ are given in Eq. (A3) . 



Finally, we note that by vector current conservation, 5Qi^ can be computed directly 
from the shifts in the up- and down-quark weak charges: 5Q^r = '25Q^^ + 5Qy^. The 
analogous relation in the SM is modified by non-perturbative strong interactions in the Z7 
box graph[6]. The latter arise because the loop contains a massless particle, rendering the 
corresponding loop integral sensitive to both low and high momentum scales. In contrast, 
the SUSY radiative corrections are dominated by large loop momenta, and non-perturbative 
QCD corrections are suppressed by (Aqcd/-^susy)^ <^ 1- 



Here, we present the SUSY one-loop diagrams that are particular to PVES. Such diagrams 
correspond to the generic corrections shown in Fig. lc,d, and e. Contributions corresponding 
to Fig. la and b are universal, and the relevant diagrams - together with the external leg 
corrections for all fermions - are given in Ref. [11]. In addition, some simplifications occur in 
the analysis for PVES that do not arise in general. In the case of charged current observables, 
for example, gluino loops can generate substantial corrections [9, 11]. In contrast, gluinos 
decouple entirely from the one-loop MSSM corrections to semi-leptonic neutral current PV 
observables. The proof of this statement is given in Appendix B. In addition, the MSSM 
Higgs contributions to vertex, external leg, and box graph corrections are negligible due to 
the small, first- and second-generation Yukawa couphngs. The fight Higgs contribution to 
gauge boson propagators has already been included via the oblique parameters, while the 
effects of other MSSM Higgs bosons are sufficiently small to be neglected [20]. Therefore, we 
do not discuss gluino and Higgs contributions in the following. 

Anapole moment corrections, corresponding to Fig. 1(c). In the presence of 
parity-violating interactions, higher-order contributions can generate the photon-fermion 
coupfing of the form (see e.g. Ref. [21]): 



where / is a fermion spinor, is the photon polarization, and F^{q^) is the anapole moment 
form-factor. 

The quantity Fj^{q^) is, in general, gauge-dependent. This dependence cancels after 



One-Loop SUSY Feynman Diagrams 




z 



(12) 
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the anapole moment contribution is combined with other one-loop corrections to the given 
scattering process [21]. The Feynman diagrams that contribute to are shown in 

Fig 2, and the analytical expressions are presented in Appendix C. 




FIG. 2: Electron anapole moment contributions to the parity-violating electron-fermion scattering 
amplitude. 

Vertex corrections, corresponding to Fig. 1(d). The relevant diagrams are shown 
in Figs. 3 and 4. The diagrams in Fig. 3 cover ee, erf, and eu scattering when the radiative 
correction is for the projectile side. When the radiative correction is to be applied to the 
target side, the diagrams in Fig. 3 can also be used for ee and ed scattering. In this case 
/ = e is the projectile. To obtain the corrections to the down quark vertex, the electron 




(a) (b) (c) (d) 

FIG. 3: MSSM radiative corrections to the electron neutral current vertex. Radiative corrections 
to the down quark vertex are obtained by replacing charged leptons (sleptons) with down type 
quarks (squarks) and sneutrinos with up type squarks. 
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(a) (b) (c) (d) 

FIG. 4: MSSM radiative corrections to the up quark neutral current vertex. 

can simply be replaced with the down quark. The diagrams in Fig. 4 show the radiative 
corrections to the target side when the incoming electron interacts with the up quark inside 
the proton. The explicit expressions for the vertex corrections can be found in Appendix C. 

Box corrections, corresponding to Fig. 1(e). These graphs generate S'q^^ in Eq. (11). 
The relevant diagrams are shown in Fig. 5 and Fig. 6. The explicit expressions are given in 
Appendix C. 



e e e e 




a) b) c) 



FIG. 5: MSSM box graphs that contribute to the electron-electron scattering amplitude. Here, 
Pi{Pi)i ^ = 1)2, is the momentum of the initial (final) state fermion. Radiative corrections to the 
electron-down quark scattering are trivially obtained by replacing the target with the down quark. 



IV. R-PARITY VIOLATING CONTRIBUTIONS TO Q 



When R-parity is not conserved, new tree-level contributions to Q^y appear. The latter 
are generated by the {B — L)-violating superpotential: 
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Wjipy — -XijkLiLjEk + X'iji^LiQjDk + -X'l-^UiDjDk + fi[LiHu , (13) 

where Li and Qi denote lepton and quark S\J{2)l doublet superfields, Ei, Ui, and Di are 
singlet superfields and the Xijk etc. are a priori unknown couplings. In order to avoid 
unacceptably large contributions to the proton decay rate, we set the AB ^ couplings 
X'-ji^ to zero. For simplicity, we also neglect the last term in Eq. (13). The purely leptonic 
terms (Ai2A;) contribute to the electron scattering amplitudes via the normalization of NC 
amplitudes to and through the definition of s^[22]. The remaining semileptonic, AL = ±1 
interactions (A^^;,) give direct contributions to the eq scattering amplitudes. The latter 
may be obtained computing the Feynman amplitudes in Fig. 7(b,c) and preforming a Fierz 
reordering. In this manner one obtains the following effective four-fermion Lagrangian: 

\y \2 _ lA' P 

^Tpv = —^^dnYdReL7,.eL + h^^L^^LeLl^^eL 

(14) 

where we have taken \q^\ <^ M| and have retained only the terms relevant for the PVES 
scattering. Note the absence from Eq. (14) of the parity- violating contact four-electron 
interaction. It is straightforward to show that the superpotential in Eq. (13) can only 
produce parity-conserving contact interactions between identical leptons. 

Contributions from P/j-violating interactions to low energy observables can be 



I \ |2 
|^12fc| 



^iiLl^lJ'LeLliil'eL + h.C. 
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a) b) c) 

FIG. 7: Tree-level Pp-violating contributions to the muon decay [plot (a)], the eu scattering am- 
plitude [plot (b)], and the ed scattering amplitude [plot (c)]. The quantities Ajjfc, eic, are defined 
in Eq. (15). 



parametrized in terms of the following quantities: 



I \ |2 



(15) 



with a similar definition for the primed quantities. In terms of Ajj^ etc. one obtains for the 
relative shifts in the weak charges [22]: 



5Q 



e 

W 



Qw 

P 
W 



1-4 sin^ e 



w 



Qw 



4 sin^ 9w 



§2(1 - s^) 



1 



1-2P 1 - AfSM 



0.35 



(16) 



As discussed in Section V the quantities Ajj^ etc. are constrained from other precision data. 
Since they are non-negative, Eq. (16) indicates that the relative shift in is negative 
semidefinite. On the other hand, the relative shift in Q^^ can have either sign depending on 
the relative magnitudes of A^k, A[-^^l., and A'^^^^.^^. 



V. ANALYSIS OF THE SUSY CONTRIBUTIONS TO THE WEAK CHARGES 

In order to evaluate the potential size of SUSY loop corrections, a set of about 3000 
different combinations of SUSY-breaking parameters was generated, chosen randomly from 
a fiat distribution in the soft SUSY mass parameters (independent for each generation) and 
in lntan/3. The former were bounded below by present collider limits and bounded above 
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by 1000 GeV, corresponding to the O(TeV) naturalness limit. Also, tan/5 was restricted 
to lie in the range 1.4 < tan/3 < 60. These limits follow from the requirement that the 
third generation quark Yukawa couplings remain perturbative (small) up to the GUT scale. 
Left-right mixing among sfermions was allowed. In order to avoid unacceptably large flavor- 
changing neutral currents, no intergenerational sfermion mixing was permitted. The ranges 
over which the soft SUSY breaking parameters and tan (3 were scanned are shown in Table I. 



Parameter 


Min Max 


tan/3 
M 

V fJLR 


1.4 60 
50 GeV 1000 GcV 
-10^ GeV^ 10^ GeV^ 



TABLE I: Ranges of SUSY parameters scanned. Here, M denotes any of ||Lt|, Mi^2, or the diagonal 
sfermion mass parameters . The ^ parameter and Mi^2 can take either sign. The generation 

index i runs from 1 to 3. 

For each combination of parameters, we computed superpartner masses and mixing an- 
gles, which wc then used as inputs for computing the radiative corrections. Only the param- 
eters generating SUSY contributions to the muon anomalous magnetic moment consistent 
with the latest results [23] were considered. We also separately evaluated the corresponding 
contributions to the oblique parameters. The latter are tightly constrained from precision 
electroweak data. We rule out any parameter combination leading to values of S and T ly- 
ing outside the present 95% confidence limit contour for these quantities. We note that this 
procedure is not entirely self-consistent, since we have not evaluated non- universal MSSM 
corrections to other precision electroweak observables before extracting oblique parameter 
constraints. As noted in Ref. [15], where MSSM corrections to Z-pole observables were 
evaluated using different models for SUSY-breaking mediation, non-universal effects can be 
as large as oblique corrections. Nevertheless, we expect our procedure to yield a reasonable 
estimate of the oblique parameter constraints. Since S and T do not dominate the low- 
energy SUSY corrections (see below), our results depend only gently on the precise allowed 
ranges for these parameters. 

In Fig. 8, we plot the shift in the weak charge of the proton, 5Qy^ — 2SQ'^ + SQ^, versus 
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the corresponding shift in the electron's weak charge, SQ^, normahzed to the respective 
SM values. The corrections in the MSSM (with Pr conserved) can be as large as ~ 4% 
(Qw) ~ iQw) ^ roughly the size of the proposed experimental errors for the two 
PVES measurements. Generally speaking, the magnitudes of SQ"^ slowly increase with 
tan (3 and decrease as SUSY mass parameters are increased. The largest effects occur when 
at least one superpartner is relatively light. An exception occurs in the presence of significant 
mass splitting between sfermions, which may lead to sizable contributions. However, such 
weak isospin-breaking effects also increase the magnitude of T, so their impact is bounded 
by oblique parameter constraints. This consideration has been implemented in arriving at 
Fig. 8. 




FIG. 8: Relative shifts in electron and proton weak charges due to SUSY effects. Dots indicate 
MSSM loop corrections for 3000 randomly-generated SUSY-breaking parameters. Interior of 
truncated elliptical region gives possible shifts due to Pr nonconserving SUSY interactions (95% 
confidence) . 

The effects of sfermion left-right mixing were studied separately. We observe that the 
presence or the absence of the mixing affects the distribution of points, but does not signif- 
icantly change the range of possible corrections. For the situation of no left-right mixing, 
the points are more strongly clustered near the origin. Thus, while corrections of the order 
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of several percent are possible in either case, large effects are more likely in the presence of 
left-right mixing. 
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FIG. 9: Various contributions to SQ^^/Q^r■. total from SUSY loops (solid line), from Sp^^^^ 
(dashed line), from 6k^^^^ (dash-dotted line), from the vertex corrections (dotted line), and from 
the box graphs (open circles). The x axis gives the lepton superpartner mass, chosen to be the same 
for both left- and right-handed first and second generation sleptons. For this graph, tan /3 = 10, the 
gaugino soft mass parameters are 2Mi = M2 = fi = 200 GeV, and masses of the third generation 
sleptons and of all squarks are 1000 GeV. For this case, SQ^r/Q^r is about half of its maximum 
possible value. The total relative correction is clearly dominated by 6k^^^^ . 



The shifts 6Ql^ are dominated by 6k,^^^^ . This feature is illustrated for in Fig. 9 
where the soft SUSY breaking parameters are chosen such that the total SUSY correction 
to Qyy is about 4%, about a half of its maximum value. For Q^r the situation is similar. 
We observe that non-universal corrections involving vertex corrections and wavefunction 
renormalization experience significant cancellations. In addition, corrections to Q'^ due to 
shifts in the ppv parameter are suppressed by 1 — 4s^. 

We find that 6k,^^^^ is nearly always negative, corresponding to a reduction in the value 
of sin2^^^(g2) = ^py{q^) sin' ew for the parity-violating electron scattering experiments 
[see Eq. (2)]. In this case, the degree of cancellation between 2J| and Qf terms in Eq. (2) 
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is reduced, yielding an increased magnitude of Q^. Since this effect is identical for both 
Q^r and Q\y, the dominant effect of 6k, produces a linear correlation between the two weak 
charges. Some scatter around this line arises from non-universal effects in (see Fig. 8). 

As illustrated in Fig. 10, within 6k^^^'^ itself, contributions from the various terms in 
Eq. (10) have comparable importance, with some degree of cancellation occurring between 
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FIG. 10: Contributions to SQ^//Q^/ from various corrections to Jk^^^'^ [see Eq. (10)]: total from 
^^SUSY line), S parameter (dotted line), T parameter (dashed line), dy^ (dash-dotted line), 

Z — J mixing and the photon self-energy (open circles), and the electron anapole moment (crosses). 
The soft SUSY parameters are the same as in Fig. 9. 



the effects of S and T. Thus, the oblique parameter approximation gives a rather poor 
description of the MSSM effects on the weak charges. In particular the quantity Sy^ in 
Eq. (10) makes a significant contribution to 6kP^^'^ . 

As evident from Fig. 8, the relative sign of the corrections to both Qy^ and Qyy - normal- 
ized to the corresponding SM values - is nearly always the same and nearly always positive. 
Since Q^y > {Q^y < 0) in the SM, SUSY loop corrections give 6Qyr > {SQ^r < 0). This 
correlation is significant, since the effects of other new physics scenarios can display different 
signatures. For example, for the general class of of theories based on Eq gauge group, with 
neutral gauge bosons having mass < 1000 GeV, the effects on Q^y and also correlate, 
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but IQ^ can have either sign in this case[6, 24]. In contrast, leptoquark interactions 
would not lead to discernible effects in Qy^ but could induce sizable shifts in [6, 24]. 

As a corollary, we also find that the relative importance of SUSY loop corrections to the 
weak charge of heavy nuclei probed with APV is suppressed. The shift in the nuclear weak 
charge is given by 5Qw{Z, N) = {2Z + N)6Q^ + {2N + Z)5Q^. Since the sign of 5Q{v/Q{y 
due to superpartner loops is nearly always the same, and since Q^r > and Q^r < in the 
SM, a strong cancellation between SQ^y and SQ^r occurs in heavy nuclei. This cancellation 
imphes that the magnitude of 5Qw{Z,N)/Qw{Z,N) is generally less than about 0.2% for 
cesium and is equally likely to have either sign. Since the presently quoted uncertainty 
for the cesium nuclear weak charge is about 0.6% [25], cesium APV docs not substantially 
constrain the SUSY parameter space. Equally as important, the present agreement of Qy^ 
with the SM prediction does not preclude significant shifts in Q'^ arising from SUSY. The 
situation is rather different, for example, in the Eq Z' scenario, where sizable shifts in Q'^ 
would also imply observable deviations of Q^^ from the SM prediction. 

The prospective "diagnostic power" of the two PVES measurements is further increased 
when one relaxes the assumption of Pr conservation. Doing so leads to the tree-level correc- 
tions to the weak charges shown in Eq. (16). The quantities Ajj^, etc. in Eqs. (15) and (16) 
are constrained from the existing precision data [22] . A summary of the existing constraints 
- including the latest theoretical inputs into the extraction of Q^^ from experiment [25] - 
is given in Table II of Ref. [11], which we partially reproduce here in Table II. We hst 
the PR-violating contribution to four relevant precision observables: superallowed nuclear 
/3-decay that constrains \ Vud\ [26], atomic PV measurements of the cesium weak charge Q^^ 
[2], the e//i ratio Re/^j. in t^i2 decays [27], and a comparison of the Fermi constant with 
the appropriate combination of a, M^, and sin^ 6yv [28]. The values of the experimental con- 
straints on those quantities are given in the last column. We also list the the PR-violating 
contributions to 5(^1 Q^y and 5Qyy/Qyy, along with proposed experimental uncertainties. 

The 95% CL region allowed by this fit in the (^Qfr/^w ^Qw/Qw plane is shown 
by the closed curve in Fig. 8. Note that the sign requirements Aijk{f), A-j^,(/) > [see 
Eq. (15)] truncate the initially elliptical curve to the shape shown in the figure. We observe 
that the prospective effects of Pr non-conservation are quite distinct from SUSY loops. 
The value of 5Qyy/Qyy is never positive in contrast to the situation for SUSY loop effects, 
whereas SQ^r/Q^y can have either sign. Note, however, that the area enclosed by the curve 
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Quantity 


A'n,(4) 






A'21fc(4) 


Value 


5\Vud?l\Vud? 


2 





-2 





-0.0029 ± 0.0014 




-4.82 


5.41 


0.05 





-0.0040 ± 0.0066 




2 








-2 


-0.0042 ± 0.0033 


SGf^/Gf^ 








1 





0.00025 ± 0.001875 




55.9 


-27.9 


-18.7 





±0.040 


SQtv/Qw 








-29.8 





±0.089 



TABLE II: PR-violating contributions to (5|Kd|VlK<il'> SQ^/Qw^ SR^/i^^ SG^/G^, SQ^^/Q^^, 
and SQ^^/Qy^^. Columns give the coefficients of the various corrections from A'-j^ and Ai2k to the 
different quantities. The last column gives the experimentally measured value of the corresponding 
quantity (for Q^, only the proposed experimental uncertainties are shown). 

corresponding to SQ^y/Q^^ > is larger than the area corresponding to SQ^,/Q^, < 0, 
implying that SQ^^/Qy^ is more likely to be positive. In addition, the magnitude of the PR- 
violating effects can be roughly twice as large as possible magnitude of SUSY loop effects 
for both Q^. Thus, a comparison of results for the two parity- violating electron scattering 
experiments could help determine whether this extension of the MSSM is to be favored over 
other new physics scenarios (see also Ref. [6]). 

VI. CONCLUSIONS 

A new generation of precise, PVES experiments are poised to probe a variety of scenarios 
for physics beyond the SM [4-6, 10]. The sensitivity of these measurements to new physics is 
enhanced because the SM values for the electron and proton weak charges are suppressed and 
because theoretical uncertainties in the SM predictions are sufficiently small[6, 10]. Here, 
we have studied the ability of these measurements to shed new light on supersymmetric 
extensions of the SM. We have observed that in a PR-conserving version of the MSSM, the 
effects of SUSY loop corrections to the electron and proton weak charges are highly correlated 
and have the same relative sign (positive) compared to the SM prediction over nearly all 
the available MSSM parameter space. This correlation arises because the corrections are 
dominated by the SUSY loop contributions to sin^ 9^-^ {q'^) - a result that would not have 
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been obvious in the absence of an explicit calculation. Moreover, the appearance of this 
correlation does not result from the adoption of any model for SUSY-breaking mediation, 
as we have undertaken a model-independent analysis in this study We also find that the 

impact of SUSY radiative corrections on the cesium weak charge are quite small, so that 
the present agreement of Qy^ with the SM does not rule out potentially observable effects 
in PVES. 

In contrast, the effects on and Q^r induced by new tree-level, Pr violating SUSY 
interactions display a different behavior. Given the constraints from other precision elec- 
troweak observables, such as the Fermi constant, first row CKM unitarity, and Qy^, one 
would expect Pr violation to cause a decrease in the size of Q^/. On the other hand, the 
magnitude of Qyy can change either way, with an increase being more likely. Moreover, the 
size of the Pr violating corrections could be even larger than those induced by SUSY loops, 
particularly in the case of Q^. Should measurements of the weak charges be consistent with 
this signature of Pr violation, they could have important implications for the nature of cold 
dark matter (it would not be supersymmetric) and the nature of neutrinos (they would be 
Majorana fermions). 

From either standpoint, should the PVES measurements deviate significantly from the 
SM predictions, one may be able to draw interesting conclusions about the character of 
SUSY. But what if both measurements turn out to be consistent with the SM? In this case, 
one would add further constraints to the possibility of Pr violation, but only marginally 
constrain the MSSM parameter space based on possible loop effects. In the latter case, the 
impact on both Qyr and is dominated by S sm"^ 9^^-^ {q'^)^^^^ . Although the projected, 
combined statistics of the two measurements would make them more sensitive to SUSY 
radiative corrections than either measurement would be independently, additional precision 
would be advantageous. In this respect, a possible future measurement of Q^y with a factor 
of two better precision than anticipated at SLAG would significantly enhance the ability of 
PVES to shed new fight on SUSYl 

* We thank D. Mack, P. Reimer, and P. Souder for sharing with us the possibility of such a future mea- 
surement at the Jefferson Lab. 
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APPENDIX A: COUNTERTERMS FOR THE EFFECTIVE PVES LAGRANGIAN 



The "bare" effective Lagrangian for the forward angle PVES scattering has the form: 



v/2 8(M0.)2 ^ ' 

Q§ = 2li - iQfsl = 27/ - iQfis' + Ss') , 

Af = (67^^756)0 = (l + ^) , 

(/V)o-^./(l + ^) ' (Al) 

where the bare quantities are indexed by "0". Unless otherwise indicated, all higher-order 
contributions include both the SM and the SUSY pieces. The quantity 1/Q{y in the paren- 
theses in the last two lines of the above equation is explicitly factored out to make the 
definitions in Eq. (A4) below more convenient. 

The counterterm 50^ is entirely determined by the muon lifetime. It can be taken from 
Eq. (62) in Ref. [17]: 



G^ 



S% , (A2) 



where flwwiQ"^) is the W boson self-energy and Sy^ is the sum of the vertex and box 
corrections to the muon decay amplitude. 
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We use the following convention for the Z-fermion interaction: 

gir = 27/-4g^s^ , 

gi = -2li . (A3) 

In this convention, the counterterms for the vector and the axial vector currents can be read 
off from Ref. [17] ^: 

5 A, = -g'y5Z% + g%5Z^ , 
SVf = g{.SZl - giSZ{ , 
5Z^^ = \ {5Zi^ - SZ^/) , 

where 5Z^£^ and 5Z^^^ are the field strength renormahzation constants for left- and right- 
handed fermions, respectively. One can write the one-loop correction to the NC vertex 
as 

-{tg/Ac)5Vf = -{tg/Ac)h, ((?(- + 75^^) / , (A5) 

where only the contributions that are not suppressed by powers of either the momentum 
transfer (-^/I^^I/Msusy) or the fermion mass (m//MsusY) are shown. The quantities G\ 
and Gy represent rescaling of the vertices by the one-loop radiative corrections. They must 
be combined with the appropriate counterterms from Eq. (A4) to obtain the renormalized 
corrections: 

VX^ + 5 A, ^G\- g-y5Z\ + g\6Z-y , 

^ Gi + SVf = Gi + g^SZl - g{6Z{ . (A6) 

APPENDIX B: DECOUPLING OF GLUINOS FROM THE WEAK CHARGE OF 

QUARKS 

It is sufficient to demonstrate the decoupling for one of the quark flavors {e.g. the up 
quark) since for other flavors the proof is identical. Consider the renormalized vector neutral 

^ Note that Ref. [17] has the opposite sign convention for g^. 
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current vertex for the up quark [see Eq. (A6)]: 



5Zl + SZ- 



5Z^ 



V 



5Zl 



5Zl - 5Zl 



(Bl) 



The gluino contributions to Gy are given by the graph shown in Fig. 4(b), with the neu- 
trahno replaced by the gluino. By using Eq. (C8) below with the appropriate coupling 
constants it is straightforward to show that: 



G^Gluino) = -\^Y. - 1^%,^ 

X '9*g19gl + 9g19gr] ^2(M^, . m^J , 



(B2) 



where as is the strong coupling constant, Mg is the gluino mass, V2(M, mi, 777.2) is defined 
in Eq. (CI) below and (see Ref. [11]) 



9gl 



9hR 



-\Pizf' 



u > 



= V2Z*^^ . 



(B3) 



In this work, no flavor mixing in the squark sector is allowed. Therefore, Z^, Z^ ^ only 
if i = 1, 4. Since Zu is unitary we find: 



9gl9gl + 9gr9gr — 2 {z]j Z^'^ + Z^ Zlf^^^ — 2Sij 



(B4) 



for i,j = 1,4. Finally: 
(Gluino) = 



_^^i9v 
3 27r\ 2 



n 

^ (1 - 2|Z^t) [V2{Mg,mij^,m^;) - V,{Mg,mi, 



774' ^Ua' 



(B5) 



where the closure property Eq. (B4) was used together with 



gl = -2/3" = -1 . 



(B6) 
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On the other hand, the gluino-induced wave function renormahzation constants of the up 
quark have the form: 



5Z^(Gluino) = l^T.i\9h'L\' + \9^'R\')Fiimu,,M^,0), 



(B7) 



5Z^(Gluino) = ^^Y:{\9h'Lf-\9h'R\')Hrriu,,M-,,Q), 
where Fi (mi, 7712,7713) is given by 

F,(™.,™„™3)= Aln{H + (l-x),„^-x(l-xK]/^^} . (B8) 

t/ 

Note that according to Eq. (CI) below, Fi(mi,m2,0) = V2(rri2, rrii, mi). Using Eqs. (B3) 
and (B4) we find 

4: as 



5Z^(G\mno) = 



3 47r 

4 as 



SZliGlvdno) = 3^ (1 - 2|Z^T) [^(M5,m^^,m^,J - y2(M5, m^^, m^j] . (B9) 

After substitution of Eqs. (B5) and (B9) into Eq. (Bl) the gluino corrections to the vec- 
tor neutral current vertex of the up quark cancel exactly. Therefore, gluino loops do not 
renormalize the weak charge of the up quark. 



APPENDIX C: COMPLETE EXPRESSIONS FOR FEYNMAN DIAGRAMS 

In this appendix we list analytical expressions for all SUSY one-loop vertex and box 
Feynman diagrams that contribute to PV electron scattering. The complete expressions for 
remaining diagrams (see Section III) are given in the appendices of Ref. [11]. We use the 
capitalized letters / and J to denote the family index for quarks and leptons {I, J — 1, • • • , 3), 
small letters i and j to denote the index for squarks and sleptons {i, j — 1, ■ • ■ ,6 except for 
sneutrino, when i,j — l,---,3), and small letters p and n to denote the index for the 
neutralinos {p,n = 1, ■ • ■ , 4) and charginos {p,n = 1, 2). 



1. Vertex Corrections 



The Feynman diagrams are shown in Figs. 3 and 4. Let us start with the corrections to 
the e — e — Z vertex. The loop integral functions Vi{mi,m2,m3) and V2(mi,m2,m3) are 
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defined as 



Vi(M, mi, 777,2) 
V2(M, 7771,7772) 

D:i{M, 7771,7772) 



/ dx dy 
Jo Jo 



D3(M, 7771, 7772) 



1 /•! r 

dx / dyy\n Ds{M,'mi,m2)/ n 

JO L 

1 — 7/)M^ + y[{l — x)ml + xm 



^1 

21 1 



(CI) 



where n is the renormahzation scale. Exphcitly: 

ml In^ 



l/i(M, 7771,7772) 

y2(-^, "71,7772) 



(M^ — rn\){ml — m\) 



21nM2 -3 + 



277i^ 



(M^ — m\){m\ — ml) 

2m\ , 777? 

in 



2\ ' 



{M^ - ml){ml - ml) 



777 



{M'^ - ml){ml - ml) 
We have [P^ = (1 - 75)72, Pi? = (1 + 75)/2]: 



(C2) 



e . 



(C3) 



Note that — and Zl^ is a unitary 3x3 matrix. Therefore, one identically has 
J2jZ*^^Zl^ — 2Qi,s^5ij — Sij. The exphcit form is kept so that the down quark neutral 
current vertex may be easily obtained by the replacement e — > d (with Zl — > Zd) and 
1/ — > 77. 



^V^^'^ = £E (r^Zl'-Zi^ + 2QjX^V,{m^o,m 
X [9oT9olPl + glTgoKPR 



(C4) 



q; 



27r^ 

+ [05:^r"^r^L + Oj;^^"^^^/^] 2777^-.777^-,yi(777~ , 777^+, 777^+) 

- o,^;^r"^r^L + o^g^rg^p^ [i + 2y2(m~ , 777^+, 777^+ 



(C5) 



+ 



Oor^?orA. + oXglTgoRp^ 2777^0777,01^(777^^,777,0,777^0) 

)R" „*ein eip f) < r^Lll *ein eip 



O^girgllPL + 0-^*r^oI^«J [1 + 2^(777^^, 777,0, 777,0) 



(C6) 
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The vector and axial vector pieces can be readily read off from the above formulae. The 
radiative corrections to the up quark neutral current vertex are as follows (see Fig. 4): 



X 



R 



U , 



_ r ^ujp uip A I *ujp uip A 

X u-f^ [9oL 9olPl + QoR 9orPr\ > 



+ 



OS9orgoTpL + Oi:;g;rgolpR] [l + 2V,{mfj^,m^o,m^o) ' 



+ 



(C7) 



(C8) 



(C9) 



(CIO) 



2. Anapole Moment Corrections 

Using formulae in the appendix of Ref. [11] we find for F|(0) of the electron: 

F1(0) = Fi")(0) + F^)(0), 
F^)(0) = 



/ . \gL 



487r 



Ff(0) 



Jo (1 - a;)m|. + xm2+ ' 



i,p 



(1 — x)m^o + xrrij 



3. Box Graphs 



Let us introduce the following notation: 



Lfi = /(p:)7m(1-75)/(a) , 
Rfi = /(pD7m(1 + 75)/(pO ■ 



(Gil) 
(C12) 



(C13) 
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As explicitly shown below, each box graph has the following structure: 

+ C'^Rf2 X Lei + D'^Lf2 X Rel] . 



(A'^L f2 X Lei + B'fRf2 x i?el 



(C14) 



To study the effects of the parity- violating electron scattering we need to pick only the term 
that has axial vector current on the projectile side el and the vector current on the target 
side /2. Therefore, it is easily seen that the box diagram contribution to Cif is [see Eq. (11)]: 



(C15) 



in the above notation. The explicit expressions for the ee box graphs in Fig. 5 are given 
below. 



SM' 



ee(a) 
Box 



^/2 47r ^^^^^^ 

I 9ol'^9ol'9ol" 901^62 X Lei + 9QR^9lR9Q'R'9mRe2 X Rel 



+ 



90R 9or9ol"' 901^62 X Lei + Qql 9TR'90RLe2 X Rel 



SM' 



ee(b) 



Box 



E 



\/2 47r ^ . 

V ^ n,p,i,3 

*ejn ejp ein *eip J j , *ejn ejp ein *eip n ^ p 

9oL 9oL 9oL 9oL ^e2 X Lei + Qqji Qqji Qqji Qqji -Ke2 X Uel 

X Bi{mi.,mi.,m^o,m^oJ 



+ 



90R 90R 90L 9oL ^e2 X ^el + 9oL 9oL 9oR 9oR ^e2 X -Kel 



X m^om^o^B2{mi^,mi^,m^o,m^oJ 



ee(c) 



Box 



—l- 



E 9r9r9r9*n2 x l, 



n,P,i,3 

X Bi{m~,m~,m^+,m^+) . 



'V2 4:71 

In all above formulae the following functions are used: 

1 rl 



el 



Bi{Mi,M2,mi,m2) = / dx [ dy f dz / 

Jo Jo Jo JJ4[Mi, 



z{l-z) 



M2,mi,m2) 



(C16) 



(C17) 



(C18) 
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B2{Mi,M2,mi,m2) = / dx dy dz -, 

Jo Jo Jo Di[Mi,M2,mi,m2) 

D^{Mi,M2,mi,m2) = z[{l - x)M^ + xM^] + {1 - z)[{yml + {1 - y)ml] . (C19) 



Explicitly: 



Si(Mi,M2,mi,m2) = 



2(m? - Mi){ml - ml){ml - M|) 

J £ 1- ^ 



2(mi-mf)(mi-Mf)(mi-M|) 2(Mf - mf)(Mf - mi)(Mi2 - M|) ' 

2 I 



B2{M„M2,m„m2) = (;;^23m^)K - -i)K " M|) 

(mi - mf ) (mi - Mf) (ml - M|) (Mf - mf ) (Mf - ml) (Mf - M|) ' ^ ' 

The box graphs for the electron-down quark scattering are easily obtained from the above 
expressions by replacing the target electron e2 with the down quark. Also, all quantities 
that have the running index j re to be replaced with the corresponding quantities for the 
first generation down squarks: g^-'j^ g^^^ etc. 

The box graphs for the electron-up quark scattering are shown in Fig. 6. The graphs (a) 
and (b) are easily obtained from the corresponding graphs for the ee scattering by replacing 
all quantities that have the running index i with the corresponding quantities for the first 
up generation squarks: g^^ — > g^Y!'-, etc. The result for the last graph (c) is: 

^ ^xj^x^^2(m~,m^.,m^+,m^+) . (C21) 



